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Title : 

CHANNEL ESTIMATION USING THE GUARD INTERVAL OF A MULTICARRIER SIGNAL 

Description 

5 Field of the invention 

This invention relates to communication using Ortliogonal Frequency Division 
Multiplexing ('OFDM') and, more particularly, to channel estimation and tracking in 
OFDM communication. 

Background of the invention 

10 OFDM communication has been chosen for most of the modern high-data 

rate communication systems (Digital Audio Broadcast - DAB, Terrestrial Digital 
Video Broadcast - DVB-T, and Broadband Radio Access Networks - BRAN such 
as HIPERLAN/2, IEEE802.11a, for example). However, in most cases the receiver 
needs an accurate estimate of the channel impulse response. 

15 In many known OFDM systems, each OFDM symbol of size JVgN* is 

preceded by a guard interval that is longer than the channel impulse response 

(CIR) and a cyclic prefix of De N+ samples is inserted as the guard interval at the 

transmitter, the prefix consisting of D samples circularly replicated from the useful 
OFDM symbol time domain samples. The cyclic prefix enables very simple 

20 equalisation at the receiver, where the cyclic prefix is discarded and each 
truncated block is processed, for example using Fourier Transform (usually Fast 
Fourier Transform), to convert the frequency-selective channel output into N 
parallel flat-faded independent sub-channel outputs, each corresponding to a 
respective sub-carrier. For equalisation purposes, numerous strategies exist. 

25 Following the zero forcing approach, for example, each sub-channel output is, 
unless it is zero, divided by the estimated channel coefficient of the corresponding 
sub-carrier. 

Like other digital communication systems, OFDM modulation encounters 
problems at high Doppler spreads, which occur notably when the user is moving 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/064344 



-2- 



PCT/EP2003/050767 



fast, for example In a car. HIPERLAN/2, for example, was designed to work only 
up to speeds of 3 m/s ("pedestrian speed"). Accordingly, the channel impulse 
response needs to be constantly tracked and updated, especially in the presence 
of high Doppler spreads. 

5 In a known OFDM communication system pilot tones are added which may 

change their position from one OFDM symbol to another. The amplitudes and 
positions of the pilot tones are known to the receiver. The receiver uses the pilot 
tones to estimate the channel coefficients of the corresponding carriers. This 
method is widely used, but it degrades the system performance, since a certain 
10 number of carriers cannot be used for data, since they are reserved for the pilot 
tones. 

It is also known to add leaming sequences (See for example EBU Review 
Technical No. 224, August 1987, "Principles of modulation and channel coding for 
digital broadcasting for mobile receiver", by M. Alard and R. Lassalle.). In 
15 HIPERLAN/2, for example, there are at least 2 leaming OFDM symbols per frame 
(i.e. 2 OFDM symbols of l Ajjs duration in total per 2ms). If the channel changes 
quickly, there must be many more training sequences and the consequence is an 
even bigger degradation in the system performance. 

Many of the known systems are unable to decode all carriers of OFDM 
20 symbols In the presence of channel nulls. Recent innovations propose ways for 
decoding OFDM symbols even In the presence of channel nulls (see for example 
the publication entitled "Reduced Complexity Equalizers for Zero-Padded OFDM 
transmissions" by B. Muquet, Marc de Coun/ille, G. B. Giannakis, Z. Wang, P. 
Duhamel in the proceedings of the Intemational Conference on Acoustics Speech 
25 and Signal Processing (MCASSP') 2000 and the publication entitled "OFDM with 
trailing zeros versus OFDM with cyclic prefix: links, comparisons and application to 
the HiperLAN/2 system" by Muquet, B.; de Courville, M.; Dunamel, P.; Giannakis, 
G, in the proceedings of the IEEE Intemational Conference on Communications, 
2000, Volume: 2. However, these publications do not offer responses to the 
30 problems referred to above concerning channel estimation and channel tracking. 

Ideally, the OFDM modulation system would keep all the advantages of 
classic OFDM and additionally allow very simple and completely blind channel 
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estimation at the receiver. No additional redundancy would be added to the 
system and therefore no bandwidth would be lost. Such a system would be 
advantageous in low-mobility scenarios and would make OFDM systems 
applicable to high-mobility scenarios as well, 

5 Many of the examples and illustrations presented below are based on the 

assumption N = 4'D, that is to say that the size of the prefix (D samples) is 
assumed to be one quarter of the size of the useful OFDM symbol {N samples). 
This con^esponds to the case of HiperLAN/2 or IEEE802.11. This restriction Is 
introduced for sake of simplicity only. It will be appreciated that the examples and 

10 illustrations are applicable more generally to the case of J\reN+ DeN* the 
necessary adaptation being basically straightfonvard. 

Summarv of the invention 

The present Invention provides a method of, and a transmitter and a receiver 
for, communication using OFDM as described in the accompanying claims. 

15 Brief description of the drawinos 

Figure 1 is a block schematic diagram of an OFDM communication system 
comprising a transmitter and a receiver in accordance with one embodiment of the 
invention, given by way of example, 

Figure 2 is a schematic diagram of an OFDM frame in a signal appearing in 
20 operation of the system of Figure 1 , 

Figure 3 is a matrix equation representing the channel impulse response for 
inter-block interference in operation of the system of Figure 1, 

Figure 4 is a matrix equation representing the channel impulse response for 
inter-symbol interference in operation of the system of Figure 1 , 

25 Figure 5 is a matrix equation representing the combined channel impulse 

response in operation of the system of Figure 1 , 

Figure 6 is a representation of a sub-matrix corresponding to the combined 
channel impulse response in operation of the system of Figure 1 for a prefix part of 
the signal of Figure 2, 
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Figure 7 represents the upper triangular sub-matrix of the channel matrix 
presented by Figure 6, 

Figure 8 represents the lower triangular sub-matrix of the channel matrix 
presented by Figure 6 , and 

5 Figure 9 is a matrix equation representing signals appearing as a result of the 

combined channel impulse response in operation of one embodiment of a system 
of the kind shown in Figure 1 , 

Figure 10 is a matrix equation representing signals appearing during channel 
estimation in operation of one embodiment of a system as shown in Figure 1, 

10 Figure 11 is a matrix equation representing signals appearing as a result of 

the combined channel impulse response in operation of another embodiment of a 
system of the kind shown in Figure 1 

Figure 12 is a graph representing preferred values of prefixes used in the 
system of Figure 1 . 

15 Detailed description of the preferred embodiments 

Figure 1 shows an OFDM communication system in accordance with one 
embodiment of the invention comprising a transmitter comprising an OFDM 
modulator 1 and a receiver comprising an OFDM demodulator 2, the transmitter 
and the receiver communicating over a communication channel 3. 

20 An input bit-stream 6„e (0,iXw = 0,1,...,a:-1 is modulated onto a set of N 

carriers whose carrier amplitudes are given by the vector 
X{k) = {X^(k),X^ikX...,Xj^_^(k)y , corresponding to OFDM symbol number k. 

Afterwards, the time domain OFDM signal is generated by means 4 which 
peri'orms an Inverse Fourier Transform operation, or preferably an Inverse Fast 

25 Fourier Transform (MFFT') operation [F^V ^U^nT with [F^Y ^iF^^j) where (J 

is the transposition operator and (•)* is the complex conjugate operator: 

X W = [F^'Xik) = ix,ikXx, Equation 1 
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Where [Fj=47 (Oi>= 



and 



N 




The resulting parallel signal x (k) vector is converted to a series signal by a 
parallel-to-series converter 5, a prefix, represented by the Dxl vector 
^i)=(^^o»" »c^.i)^. being inserted into the signal as guard interval between each 
5 OFDM symbol to produce a series digital signal x„ . The series digital signal x„ is 
then converted to an analogue signal x(0 by a digital-to-analogue converter 6 and 
transmitted over the channel 3. 

The channel 3 has a Channel Impulse Response H(k) = C(k) and also 
introduces noise v. 

10 At the receiver 2, an analogue signal r(t) is received and converted to a 

digital signal r„ by an analogue-to-digital converter 7. The digital signal r„ is then 
converted to a parallel signal by a series-to-parallel converter r(Jfc) and equalised 
and demodulated by equalisation and demodulation means 9 to produce 
demodulated signals s^'(k) . In the following analysis, consideration of noise is 

15 omitted for the sake for simplicity. However, including the consideration of noise 
does not significantly modify the results. 

In some known OFDM communication systems, the guard interval is used to 
add some redundancy {D samples of redundancy are added) by introducing a 
cyclic prefix, for example in the following manner: 



In other words, data from the end of the frame is repeated by the transmitter in the 
guard interval to produce a prefix. 

In accordance with this embodiment of the present invention, however, the 
prefix samples inserted as guard interval of OFDM symbol number k, Ok.co to a^.co. 
25 t, are deterministic and are known to said receiver as well as to said transmitter. 
The prefixes comprise a vector i'o =(^ov M^z).i)^of size Z)xl that is common to the 

symbols multiplied by at least one weighting factor aj^, so that the prefixes have the 
overall form qi^.cq to ai^.CD-i- The weighting factor may be constant from one 



20 x''^'{k) = ix,^M 
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symbol to another. However, in a preferred embodiment of the invention, the 
weighting factor Ok differs from one symbol to another, the elements of a given 
vector Pq being multiplied by the same weighting factor. With an OFDM modulator 
in the transmitter functioning in this way, blind channel estimation in the receiver 
5 can be done simply and at low arithmetical complexity. In particular, the receiver 
can constantly estimate and track the channel impulse response without any loss 
of data bandwidth. Moreover, the demodulator at the receiver can have 
advantageous characteristics, ranging from very low arithmetical cost (at medium 
performance) to high arithmetical cost (very good system performance). 

10 More partlculariy, in the preferred embodiment of the invention, the prefix of 

D samples that Is added in the guard interval comprises a pre-calculated suitable 
vector i'/) =(Co,...,c^.,)^ of D samples that is independent of the data and that is 
weighted by a pseudo-random factor a^^ that only depends on the number k of the 
latest OFDM symbol: 

15 x^^'\k) = (a,Co,...,a,c^.pXo(^),:c, {kU,x^Jk)y, Equation 2 

For the purposes of the analysis below, a second prefix/OFDM symbol vector 
is defined as follows: 

^ico....poso^f^^ = ixM^ (4....%-,(4^*^i^ov..,a*+iC^.,/. Equation 3 

Several choices for are possible. It is possible to choose a,^ e C, that is to 

20 say that can be of any complex value. However, any a^wlth \aj^\^l leads to 
performance degradation compared to preferred embodiments of the invention. 

It is possible to limit the choice of , somewhat less generally to a^^C with 

= This choice usually leads to good system performance, but the decoding 
process risks to be unnecessarily complex. 
25 Accordingly, in the preferred embodiment of the present invention, the phase 

J m 

of aj^ is chosen so that o:f^=e , where m is an integer, N is the useful OFDM 
symbol size and D is the size of the pseudo-random prefix. This choice is 
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particularly advantageous when using the specific decoding methods described 
below. 

For the sake of simplicity, the following analysis assumes that the weighting 

J m 

factor has been chosen as aj^ = e ,m integer. However, it will be appreciated 

5 that the mathematical adaptation to any of the cases presented above is 
straightforward. 

It proves to be very useful to choose «^ such that its phase changes from 
OFDM symbol to OFDM symbol. The constant prefix Pq is preferably chosen with 
respect to certain criteria, for example the following: 
10 • In the frequency domain, is as flat as possible over the frequency 

band used for data carriers. 

• In the frequency domain, is as near to zero as possible for all 
unused parts of the band. 

• In the time domain, P^, has a low peak-to-average-power-ratio 
15 (PAPR). 

• The length of P^ is the size of the OFDM guard interval, that is to say 
D samples. Alternatively, a shorter sequence of length b<D may be 
chosen where D-b zeros are appended. 

With these criteria, without any complication of the transmitter, the receiver is 
20 able to estimate the channel impulse response blindly, track the changes of the 
channel impulse response blindly and perform an arithmetically simple 
equalization. 

An example of a frame of OFDM symbols in accordance with a preferred 
embodiment of the invention is illustrated in Figure 2. The operation of the system 
25 will first be described for the specific case where «^ is constant and equal to 1 . 

Now, the modulation unit of the transmitter Is clearly defined. In the following, 
the operations to be performed in the receiver are considered. Each received 
OFDM symbol selected at the input of the demodulator 9 can then be expressed 
as follows (neglecting additive noise): 
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Kk)^[H,J^x^'^'^^^'^{k^l)+[H^,]x^^^^^ Equation 4 

where the channel impulse response of the demodulator 9 is assumed to be 
/i = (/io, [ffiBf] "s the contribution of the demodulator 9 channel matrix 

corresponding to inter-block-interference and [h^] is its contribution to Inter- 
5 symbol-Interference. 

The components of the received signal corresponding to inter-block- 
interference [^/5/lr(Ar-l) are illustrated in Figure 3, where blank elements 
correspond to zero values, for an example where N^A D (for example, m the 
case of HiperLAN/2 or IEEE802,11, iV = 64 and D = 16). It will be seen that [Hj^j] 
10 is a matrix of size {N'^D)x{N^D) with a triangular sub-matrix [h^] of size (D-^)x(D- 
1) at its upper right-hand comer, illustrated by Figure 7, the other elements of the 
matrix being zero. 

The components of the received signal corresponding to inter-symbol- 
Interference are illustrated in Figure 4, for the same case and In the 

15 same manner as Figure 3. It will be seen that [h^] is a matrix of size 
(N+D)x(A/+D) with triangular sub-matrices [if J on its major diagonal as illustrated 
by Figure 7 and triangular sub-matrices [i/J of size DxD on the diagonal 

Immediately below the main diagonal as illustrated by Figure 8, the other elements 
of the matrix being zero. 

20 The channel impulse response seen by demodulator 9 is represented by the 

sum of the inter-block-interference [if^^J and the inter-symbol-interference [Hjsj], 
as shown in Figure 5. The resulting signal for this example is shown in Figure 9, 
where r^^k) to r^(k) are successive parts of the OFDM symbol #k containing as 
well contributions of the preceding and following prefix convolved by the channel, 

25 x^ik) to Xy(k) are corresponding parts of size D of the useful signal transmitted 
and x^ik) is a corresponding part of size D of the following prefix in this example. 
Of course, the example may be generalised to any e N+, /) e N+. 

The expectation values of the parts of the received signals are as follows: 
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Eo = Eir,)=E([H,] xM+E([H,] x^{k)) = [H,] p^ Equation 5 



=0 



E, = E{r,)= E([H,] x,{k))+E([H,]'X,{k))^0 Equation 6 



^2 = Eir^) = E([H^]•x^{k))+ e([HoI x^{k)) = 0 Equation 7 
> , ' < , ' 

=0 =0 

£3 = Eir^) = E([Ht]- Xi{k))+E([Ho] • x, {k)) = 0 Equation 8 



E, = E{r,) = E([H,lx,{k))+E([H,] x,{k)) = [Hj PD Equation 9 

^ V * 

=0 

It will be appreciated that the expectation values of the useful parts x^^k) to 
x^{k) of the OFDM symbol tend to zero over a large number of symbols since they 
are quasi-random with zero mean. However, the prefix is known to the receiver 
(and in this embodiment is constant over successive symbols) and enables 
[h J= [iyo]+[^i] to be estimated, by approximating the expectation values and 
over a large number R of symbols: 

E{r) « ^ • XKO- Equation 10 
The sum of the expectation values E^ and E^ is then given by: 

E,+E,= ([H,h[H,])'P^^[HlP^ Equation 11 

A first embodiment of a method of channel impulse response estimation on D 
symbols in accordance with the present invention utilises the expression of the 
above equation as follows: 



Equation 12 
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Where the matrices [fJ and [/^^^ =[F^f = ([F^r)* are the (Fast) Fourier 
Transform and Inverse (Fast) Fourier Transform matrices respectively and the 
prefix is of size D. The matrices [i/ol[J^i] and [h] are Illustrated by Figure 7, 
Figure 8 and Figure 6 respectively. 

Accordingly, in this first method, the channel impulse response is estimated 
using the following steps: 

• Perform a FFTo^ on V^^ = [h,])- P^^E^-^ 

• Perform a FFT^^ on = P^ 

• Perform a component-by-component division of the first result by the 
second H^^^ = Vj^p 0 F> : 

fl^ir> = Fi^r^,J([iyJ+[//,]).p^)® FFT^M- 

• Perform an IFFT on H^J"^ : 



The resulting channel estimation is hj^ of size Dxl. This method works well 
in many circumstances and has a low arithmetic cost, since its calculations are 
based on matrices of size DxD. However, an OFDM symbol which usually is of 
size N>D samples will be equalized based on this estimation. Thus, this method 
worlds very well if the prefix-spectrum is non-zero everywhere in the FFT^^ 
domain (and, of course, everywhere well above channel noise). This can be a 
troublesome limitation in other circumstances. 

A second embodiment of a method of channel impulse response estimation 
on D carriers in accordance with the present Invention avoids this limitation, at the 

expense of increased arithmetic cost. This second method does not estimate 
based on a de-convolution In the FFT^ domain as presented above, but 



estimates 




directly based on the received vector 



([ifo] + [^i])'^o • This is possible by exploiting the observation: 



(A^+J5)x(iV+Z>) 



]• (PJO^^) = (e/eJ01_J Equation 13 
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This equation is represented in more detail in Figure 10. in this second 
method, the channel impulse response is estimated using the following steps: 

• Perfomi a FFT^^^j^^^o) on ^hf = [^(a^+dm^+o)]- (P/^lJ 

• Perform a FFT^^^om/v+d) on V, = (pjol J 

5 • Perfomn a component-by-component division H^/^^ = V^p ® Vp 

• If desired, perform an IFFT on H^P^^^ : ^^^^^ = IFirr^^^^y^^^^^^(H^^^^^^ 

The last step of the list presented above is not essential for the basic 
equalization algorithm but may be useful, for example in algorithms used to reduce 
noise levels. 

10 The above methods have been described with reference to the specific case 

where is constant and equal to 1. In preferred embodiments of the invention, 
however, the weight of the prefix to each symbol A: is a preferably complex 
pseudo-random factor that only depends on the number k of the latest OFDM 
symbol. The adaptations to this method of the basic equations (shown in Figure 9) 

15 are shown in Figure 1 1 . 

It is found that equations 4 and 8 are to be adapted as follows: 



=0 



E,^,=E([H,]r,{k))+E([H,l^^^^^ Equation 15 

' 3 ' V ^*+! J 

The procedures for blind channel estimation described above remain 
20 applicable by setting Eq-E^q and E^=^E^^,. This amounts to weighting the 
preceding and following D prefix-samples of each received symbol by the 
corresponding a^^ or aj^li respectively. 

The values of the prefixes aj^ Po are chosen as a function of selected 

criteria, as mentioned above. Values that have been found to give good results 
25 with the criteria: 
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• Low Peak-to-Average-Power-Ratio of the time domain signal 

• Low Out-of-Band Radiation, that is to say maximise the energy of the prefix 
over the useful band and not waste prefix energy over null carriers 

• Spectral Flatness, e.g. SNR of each channel estimates shall be approx. 
5 constant 

• LoW'Complexity Channel Estimation, i.e. by prefix spectrum whose spectral 
contributions are mainly just phases (i.e. of constant modulus), 

are shown in Figure 12 by way of example, for the following OFDM parameters: 

• Size of the Prefix in Time Domain: D=16 Samples 

10 • Size ofthe OFDM symbols in the frame: A^64 Samples 

• Carriers where channel coefficients are to be estimated (over A/+D=80 
carriers): Carriers 1 to 52 

• Out-of-Band region: Carriers 76 to 80 

• Maximum PAPR has not been limited 

15 • Out-of-Band Radiation as low as possible 

• Spectral Flatness as good as possible. 

The channel estimation is done by calculating the expectation value over a 
number of samples of the received vector as explained above. If the tracking of the 
channel is done based on a first estimation h(k - 1) = [f^]^ • H{k - 1) of the channel 
20 impulse response and a number R of OFDM symbols, the first estimate is then 
updated as follows: 

based on the ideas of the flrst method for channel estimation that has been 
25 presented above. Alternatively, the second method can be applied by 

H(k) = s, • Hik - 1) + (^[iv. J- Z-^„.. • (Elik - ri),Elik - «),0L j) O \f,,,]- {Pl,Oll] 
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where the factors 5„,n = 0,l,...,/?-l are positive real numbers that are used for 
normalization and weighting of the different contributions. Thus, for example it is 
possible to take older OFDM symbols less into account for the channel estimation 
than later ones. The Fourier matrix [f] can be chosen in the A/+D carriers or D 
5 carriers domain 

Several equalization methods are advantageous using the pseudo-random 
prefix OFDM. In general, the different methods offer different performance- 
complexity trade-offs. 

A first embodiment of a method of equalization uses zero forcing in the A/+D 
10 Domain and offers low complexity equalization. 

With >9jt = , the Channel Impulse Response matrix can be represented 
as follows: 



[h] =[h,J+-^^ [h,J = A • [h,J 



a. 



Equation 16 



Still assuming that the length of the channel impulse response is D, the 
15 coefficients /id+p -»^Ar+z)-i ^^e set to zero. This is a so-called pseudo-circulant 
matrix, corresponding to the case where is not equal to 1, and can be 
diagonalized as follows: 



JJ 



20 



1 

Where M = l S |A| ^-^^1 



and H(z)= J^z-" ■}!„ Equation 17 



n=0 
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10 



15 



With the assumption that r{k) = [Hl{x{kY P^J , and that the weighting factor 



has been chosen as A = = ' 



, where m Is an integer, and if the received 



vector R{k) is: 

the procedure of this method of zero forcing equalisation Is: 
• Perform Multiplication R^^\k) = 4nTd • [F] • R{k). , 



r 1 f ^ jv+D-i^l 



Where [f]= £ lAl"!^ 

• Calculate the frequency shifted, estimated CIR coefficients 



• Perform a component-by-component division R^^\k) = Q JHf^f^ 

• Perfonn ly/lultiplication J?(3>(it) = -,=i= — [vY-R^^\k). 

• Extract tlie N equalized samples of tlie /cth OFDM-data symbol to S^Q(/f). 

• Transform tfie kth OFDM data symbol S£Q(/f) into the frequency domain by 
a Fourier transform S^Q(/f)=[Fw] S£Q(/c). 

• Proceed witti metric calculation, etc. on tlie received equalised carriers. 

Anotlier embodiment of a method of equalization uses a method known from 
studies on zero padding. The received vector R(k) in the OFDM Pseudo Random 
Prefix Scheme can be expressed as follows, where [P\ contains a {N+D)xN pre- 
ceding matrix and In is the NxN identity matrix: 



20 



R{k) = [H]- [P] x{k)+ 



= [Hi 



Equation 18 



The Channel Impulse Response estimation obtained as described above is 
then used together with the known values of Po to perfonn the following operation 
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10 



15 



Equation 19 



in which the known prefix values are multiplied by the Channel impulse Response 
estimation and the result subtracted from the received signal. In the general case, 
[H\ is a pseudo circulant channel matrix. So, the diagonalisation of such matrices 
can then be performed in order to calculate efficiently. Then, several 

equalization approaches are possible, for example the Zero Forcing (ZF) approach 
or the Minimum Mean Square Error (MMSE) equalization approach. Examples of 
MMSE equalization methods are described in the articles "OFDM with trailing 
zeros versus OFDM with cyclic prefix: links, comparisons and application to the 
HiperLAN/2 system" by Muquet, B.; de Courville. M.; Dunamel. P.; Giannakis, G. 
ICC 2000 - IEEE International Conference on Communications, Volume 2, 2000 
and "Reduced Complexity Equalizers for Zero-Padded OFDM transmissions" by 
Muquet, B.; de Courville, M.; Giannakis, G. B.; Wang, Z.; Duhamel, P. 
International Conference on Acoustics Speech and Signal Processing (ICASSP) 
2000. 

In one example, the equalisation is performed based on a zero-forcing 
approach by multiplying y^^^ by the Moore-Penrose pseudo-inverse [g] of the 



matrix [h]' 



[Hi 



Thus, the equalized resulting vector is 



= [G]- 



[G]R'Hk) 



Equation 20 



20 The definition of tiie Moore-Penrose pseudo-inverse is, among others, 

discussed by Hayl<in in the book: ''Adaptive Filter Theory^' by Simon Haykin, 3'"^ 
edition. Prentice Hall Information and System Science Series, 1996. Haykin uses 
the common definition 



[aY = (a"aYa". Equation 21 



SUBSTITUTE SHEET (RULE 26) 



wo 2004/064344 PCT/EP2003/0S0767 

-16- 



where [a] is a rectangular matrix. 
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